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Abstract 

The generalized pseudospectral method is employed to study the bound-state spectra of some 
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of the exponentially screened Coulomb potentials, viz., the exponential cosine screened Coulomb 



(ECSC) and general exponential screened Coulomb (GESC) potential, with special emphasis on 
higher states and stronger interaction. Eigenvalues accurate up to eleven significant figures are 



obtained through a non-uniform optimal spatial discretization of the radial Schrodinger equation. 



All the 55 eigenstates of ECSC potential with n < 10 and 36 eigenstates of GESC potential with 
n < 8 are considered for arbitrary values of the screening parameter, covering a wide range of 
interaction. Excited states as high as up to n = 18 have been computed with high accuracy for 
the first time. Excellent agreement with the literature data has been observed in all cases. All the 
GESC eigenstates are calculated with much greater accuracy than the existing methods available 
in literature. Many / ^ states of this potential are reported here. In both cases, a detailed 



variation of energies with respect to the parameters in potential is monitored. 
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I. INTRODUCTION 

Realistic potentials which describe quantum mechanical systems, are not usually exactly 
solvable in the Schrodinger picture, except for a few occasions such as Harmonic oscillator, 
Coulomb potential, etc. Therefore, finding exact analytical solution of Schrodinger equation 
for a given potential corresponding to a physical system of interest, constitutes one of the 
major challenges in quantum mechanics. This is a common problem, and often encountered 
in almost every branches, such as atomic, molecular, solid-state, nuclear, particle and plasma 
physics, etc. A large number of attractive promising approximate formalisms have been 
developed ever since the inception of theory, which can provide highly accurate or near- 
exact results in some cases. However, the same for a general potential for any allowed 
values of quantum numbers, for arbitrary values of potential parameters (if present in the 
system) still remains elusive, and thus always has been an active area of research. 

Here we are concerned with the accurate bound states of two central singular potentials, 
namely (i) the generalized exponential cosine screened Coulomb potential, given by, 

A A 
v(r) = e" lT cos(g52r) = e~ r cos(g6r) (when S± — 82 — 5), (1) 

and (ii) a much less frequently studied, general exponential screened Coulomb (GESC) 
potential of the form, 

v(r) = -- [l + (l + 6r)e- 26r ]- (2) 

In Eq. (1), A represents the coupling strength constant while 61,62 are two screening pa- 
rameters. This potential reduces to the familiar Yukawa potential for g — 0, which has wide 
applications in nuclear, solid-state and plasma physics. For g — 1, this is termed as the expo- 
nential cosine screened Coulomb (ECSC) potential, and it is in this form that this potential 
has been studied maximum. One distinctive feature of this oscillating potential (in contrast 
to the Coulomb potential) is the finite number of bound states, i.e., such states exist only 
for certain values of the screening parameter below a threshold limit (the so-called critical 
6 C ). In other words, the total number of different energy levels is finite for a given value of 
6 > 0. Similarly, the two potential parameters a, b in Eq. (2) signify coupling strength and 
screening parameters respectively. Throughout the whole article, a is fixed at unity. These 
two potentials can be used to represent the effective interaction in many-electron atoms; also 
they have important applications in solid-state, nuclear and plasma physics as well as in field 



theory l|-|9|. Lately, the effect of screening on atomic pho toionization in H and He + has 
been studied by means of Yukawa and ECSC potential 10] . Also, the ground and excited 
resonances in two-electron systems such as He, molecular H 2 in ECSC as well as generalized 
screened potential have been investigated llMl4|. 

None of these potentials admits exact analytical result. Therefore, over five decades, 
a considerably large number of attempts have been made to calculate their eigen spectra 
accurately Here we mention a few of them. Perturbation and variational methods were 
used [151 to produce eigenvalues with reasonable accuracy, as well as the number of bound 

nn 

states for a given value of 5. The s states were reported [16J, LL7J via the representation 
of ECSC potential by a Hulthen potential with an energy- dependent strength parameter, 
through the use of an Ecker-Weizel approximation. The n < 4 states have been calculated 
by means of a hyper- virial Pade approximation 18] . dynamical group approach [l9j, hyper- 



20|, etc. Using a numerical method 



2l|, 



virial equation with Hellman-Feynman theorem 
all the 36 states below n < 8, as well as the critical screening parameters were obtained 
within an accuracy of eight to six significant figures. Analytical expressions for eigenvalues 
and eigenfunctions of ground and first excited states, up to 14 terms, were presented by a 
large- N expansion method [22]. Eigenvalues of Is to 8k states have been reported by means 
of a shifted 1/N expansion technique [23] . An iterative solution for eigenvalues belonging 
to arbitrary n, I quantum numbers has been put forth by employing an asymptotic iteration 
method 24J. A novel perturbation method 25] has also been proposed for this potential 
where the radial Schrodinger equation is decomposed into two parts, one of them being 
exactly solvable while the other part leading to closed analytical solution or an approximate 
treatment depending on the potential in question. Lately, a J-matrix approach 26] with 
a Gaussian quadrature scheme has offered high-quality results for bound and continuum 



states. A Ritz variation method with hydrogenic wave function as the trial function 27 1 
has also produced promising results for such potential. Recently, an analytical scheme 28] 
inspired by the J-matrix method, has been quite successful for such potentials. Th e S r values 



have been estimated by numerical 
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261 ] . as well as analytical methods 29 1. 



The GESC potential, on the other hand, has not received much attention. I am aware of 
only two studies. Energy eigenvalues of ground and first excited states were presented up to 
1 4 terms using a large- N expansion method [30] . In another attempt 3l[, the perturbative 



method of 



25] was used to obtain the Is, 2s, 3s states of this potential with decent accuracy. 



In this work, we study the eigenspectra of both these potentials in Eqs. (1) and (2) using 
a generalized pseudospectral (GPS) method, which has been quite successful for a variety of 
systems such as the spiked harmonic oscillator, Hulthen and Yukawa potentials, power-law 
and logarithmic potentials, ground and excited states (low- and high-lying Rydberg states) 



of atoms as well as other singular systems 32M36I. l38l. l39l] . Potential parameters are scanned 



over a large domain. In few occasions, for some of the methodologies mentioned above 
for ECSC potential, it so happens that, eigenvalues and eigenf unctions are quite difficult to 
calculate for high-lying states and also at certain region of the potential parameter, especially 
near the 5 C . Here we pay special attention to both these issues for a better understanding of 
their spectra and also to judge the validity and efficacy of the method. To this end, accurate 
energies and wave functions are presented for all the 55 levels belonging to n — 10 states of 
ECSC potential to extend the domain of applicability of the GPS procedure. Variation of 
the same with respect to 5 are also monitored. For the GESC potential, only some low-lying 
I = states have been considered so far in the literature; no results are available for I ^ 
states. So here, we report all the states up to n = 8 for a wide range of parameters in the 
potential for the first time, and some higher states as well. A detailed comparison with the 
available results in the literature has been made, wherever possible. The article is organized 
as follows. An outline of the theory and method of calculation is presented in Section II. A 
discussion of the results is given in Section III, while a few concluding remarks are made in 
Section IV. 

II. THE GPS METHOD 

The section gives the essential steps of GPS approach, as implemented here for the 
solution of single-particle Schrodinger equation for a non-relativistic Hamiltonian containing 
an exponentially screened potential term. The key advantage of the approach is that it offers 
a non-uniform, optimal spatial discretization. That means one can use a finer grid at small 
r and coarser grid at large r, maintaining high- accuracy at both these regions. This also 
implies only a small number of spatial points suffices to achieve convergence. Thus compared 
to standard finite difference/finite element methods, the GPS scheme is both accurate and 



efficient. Other details could be found in the references 
atomic unit is used throughout the article. 
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. Unless otherwise mentioned, 



The radial Schrodinger equation can be written in the following working form, 

1 d 2 £{£+!) 



H h v(r) 

2 dr 2 2r 2 l j 



#n,*(r) = E n/ R ny e{r) (3) 



where v(r) is as given in Eq. (1) or (2), whereas n, £ signify the usual radial and angular 
momentum quantum numbers respectively. 

As a key step, a function f(x), defined in the interval x G [—1,1], is approximated by the 
N-th order polynomial, /jv(x), through a cardinal function gj(x), as follows: 

N 

f(x)^Mx) = '£f(x j )g j (x). (4) 

3=0 

This guarantees that the approximation is exact at the collocation points Xj, i.e., Jn{xj) = 
f(xj), and requires that the cardinal function satisfies gj(xj>) = Sj/j. Here we use the 
Legendre pseudospectral method, where xq = — 1, x^ = 1, and the Xj(j = 1, . . . , N — 1) are 
obtained from roots of first derivatives of the Legendre polynomial, Pn{x) with respect to 
x, as P' n (xj) = 0. The gj(x) are given by, 

a . (x) l a-* 2 ) no*) (5) 

9A) N(N + l)P N ( Xj ) x-Xj ' {> 

Now, the semi-infinite domain r£ [0, oo] is mapped onto a finite domain x G [—1,1] via the 

transformation r = r(x). At this stage, one could introduce an algebraic nonlinear mapping 

of the form, 

r = r(x) = L ■ — , (6) 

1 — x + a 

with L and a = 2L/r max as two mapping parameters, to obtain a transformed differen- 



tial equation as: f(x) = R n i(r(x))/Jr'(x). Now, one applies the Legendre pseudospectral 
method to this equation and finally a symmetrization procedure to yield the following sym- 
metric eigenvalue problem, 

iV-l r 

~ 2 



H -~ D iJ+ u i 6 iJ Xj = tmXi- (7) 



This is readily solved by standard available routines such as that in NAG Fortran library, 
giving highly accurate eigenvalues and eigenfunctions. After some straightforward algebra, 
one finds that, 

Xi = R n i(r i ) y/W)/P N (xi), u t = l(l + l)/2^ + v(r t ), (8) 



with Xi = x( x i)i u i = u ( x i), r i = r ( x i), r i = r '( x i), while Dy denotes the symmetrized second 
derivative of cardinal function, given as follows, 

2 

V = 77 \o / ' l I 3i 

< iV l 3 ) j 



N(N + 1) 



J- (9) 



3f?(l-s?)' 

A series of calculation was performed for various potential parameters with respect to the 
grid mapping parameters to ascertain the accuracy and reliability of the current method. In 
this way, a "stable" grid was found, which appears to be sufficient for all the converged results 
presented in this article. Unless otherwise mentioned, all the reported results correspond 
to this consistent set of parameters, a = 25, N = 200 and r max = 300. There are some 
instances, where this set is not adequate, and appropriate variations are allowed; these are 
mentioned appropriately in the text. Current results are reported only up to the precision 
that maintained stability. Eigenvalues are truncated rather than rounded- off, and hence may 
be considered as correct up to all the decimal places they are reported. 

At this stage, a few remarks may be made regarding the GPS method. Typically in direct 
numerical methods, one truncates the semi-infinite domain into a finite domain [r min , r max ] to 
deal with the problems of singularity at r = 0, and infinite domain. In order for this, r min and 
r max need to be chosen sufficiently small and large respectively. This consequently results in 
a rather large number of grid points and also, in general, introduces some truncation error. 
To overcome this problem, one can map the semi-infinite domain [0, oo] exactly into the finite 
domain [—1,1] using the mapping r = f(x) (Eq. 6) so that the Legendre pseudopotential 
technique can be applied. This introduces an additional undesirable feature; namely it leads 
to an unsymmetric or generalized eigenvalue problem, which in turn, is bypassed via the 
symmetrization procedure mentioned above. The method has been successfully applied to 
resonance states as well. For these and many other features of the method, the interested 



reader is referred to the references 
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421 ] and those therein. 



III. RESULTS AND DISCUSSION 

At first, in Table I, we report some s states of the ECSC potential for low as well as high 
excitations. For all the ECSC potential calculations throughout the article, parameter A is 
set to unity. A wide range of screening parameters is considered-including low, intermediate 



TABLE I: Calculated negative eigenvalues E (in a.u.) of some selected s states of the ECSC 
potential for various 5 along with the literature data. The 5 C values of Is, 2s, 3s, 4s, 8s are 0.72, 
0.16656, 0.0724, 0.0404 and 0.01 respectively, taken from [26]. 



State 



—Energy 



-Energy 



This work 



Reference 



This work 



Reference 



Is 0.0002 0.49980000000 0.499800° 



0.1 



c , 0.499899" 



0.40088477464 0.400883 a ,0. 400884°, 0.400885°'°'™, 
0.402155 d ,0.40088477 / , 0.4008839 9 '' 1 , 
0.40088421 ! ,0.4008847&', 
0.4008785*,0.400884774639' 

2s 0.06 0.06742110520 0.067385 a ,0.067408\ 0.067421°-°, 

0.067525 d , 0.0674217° - h , 0.06742608 s , 
0.06742085 J ,0.0673900 fc , 
0.06742110514 i ,0. 06742173" 

3s 0.04 0.01882306336 0.018707°,0. 018768°, 0.018822°, 

0.019604 d ,0.018823 e , 0.0188478 s , 
0.01886716\0.01882rJ, 
0.0188586*,0.01882306333' 

4s 0.0005 0.03075002676 0.030750°'°'°, 0.030751 d 

4s 0.02 0.01257177727 0.012539°,0. 012557°, 0.012572°'°'^, 

0.013084 d ,0.012581is, 0.01259233 1 

8s 0.0001 0.00771250340 0.007713 m 
17s 0.0005 0.00123778635 
18s 0.0005 0.00105272559 



0.06 



0.7 



0.44020051029 0.440200"'' 1 '*, 0.440201°'°'°, 
0.44020051^, 0.44020049, 
0.44020057 B ,0.44020051029' 

0.00115044274 -0.050624°, -0.036908°, 
0.00184°, -0.000043°, 



0.00115044272 



0.165 0.00018502068 0.0001850' 



0.072 0.00009790825 0.0000979' 



0.005 0.02627512430 0.026275° 



°, 0.026321° 



0.04 0.00014026953 -0. 001079°, -0.000670°, 
0.000118°,0.000125°, 
0.00106945,0.00032273* 

0.005 0.00314139349 0.003134™ 

0.001 0.00078491417 

0.001 0.00061036597 



"■Perturbation (Coulomb), Re f. IE 
d Ref.Jl7|]. °Ref. JlS 

j Ref. bi. k Ref. [25 



b Perturbation (Hulthen), Ref. 11511 - c One-parameter variational, Ref^ 115 

f Ref. J21J. g Ref. [19j|. h Ref. [22j. s Ref. i23 

'Ref. [26J. m Two-parameter variational, Ref. tl5j. n Ref. [2, 




and high values, signifying small, intermediate and large interaction respectively. Critical 
values of the screening parameter, taken from [26j , are also mentioned in the table for n < 8 
for easy understanding. As evident, a large number of results are available in the literature 
for comparison, which we quote accordingly. One of the very first definitive calculations 
of this potential was reported in 15| . All the s states considered here with n=l,4 (except 
2s, 3s in the high-screening region) were estimated by first-order perturbation treatment 
with (a) Coulomb potential as unperturbed potential (b) Hulthen potential as unperturbed 



potential, and (c) a one-parameter variational calculation with reasonably good accuracy. 
An Ecker-Weizel approach has been used for ls-4s states in the low-screening region through 
an approximation of the ECSC potential by Hulthen potential with modest accuracy 17 ]. 
Both s and I 7^ states of n < 
hyper- virial Pade approximation 



were treated in intermediate coupling; region within the 



iteration method 



18], a dynamical group approach [19], and an asymptotic 



241 ] . The ground and first excited states of the ECSC potential for medium 



values of screening parameter have been obtained within the large- N expansion method 22] 
as well. Energies correct up to six to eight significant figures were reported by means of a 
shifted 1/N expansion 23] , for both I = as well as I ^ states. Lately, a new perturbative 
scheme [25] has been put forth for the n < 3 states with decent success. However, it seems 
that, so far the most accurate eigenvalues are reported by a J-matrix method [26]. In the 
neighborhood of low and moderate coupling, the present results are of very similar accuracy 
as those from 26j (in many occasions they coincide; otherwise they differ in the 11th or 



12th place of decimal). In some states, for <5s near the threshold limit, their results were 
reported for somewhat lesser accuracy. Through the present method, we are able to obtain 
eigenvalues of consistently better accuracy near the strong- coupling region (see, for example, 
2s and 3s for 5 = 0.165, 0.072 respectively). For higher-lying states, the reference literature 
values dramatically reduce in number so much so that for n = 8, our results could only 
be compared with the lone two-parameter variational calculation 15], where the present 
results are visibly improved. And for states with n > 8, no results could be found, and we 
report some sample results for 17s and 18s to emphasize the ease and ability of our method 
for higher states. Note that, for these higher states, however, the r max value needs to be 
suitably increased to achieve the desired convergence, e.g., for 17s, 18s, an r max = 1100 
a.u., was used, whereas the other two parameters a and N needed no adjustments. This is 
reminiscent of a situation encountered earlier for Hulthen and Yukawa potentials 35] in the 
stronger coupling regime. Some numerical results are also available for the low-lying states 
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27] . which have been quoted as well. 



Next in Table II, we report all the / 7^ states belonging to n = 2 — 6, at selected values 
of S. Here also the 5 values are chosen so as to reflect both weak and strong couplings of 
the interaction, with respective 5 C values mentioned at the top of table. While the literature 
results are clearly quite scanty in comparison to / = case in Table I, wherever available, 
these are quoted appropriately. Once again, the existing best result is apparently the one 
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TABLE II: Calculated negative eigenvalues (a.u.) of the ECSC potential for selected n = 2— 6,£ 7^ 
states for various 5 values along with the literature data. The S c values for p, d, f, g, h states for 
n = (I + 1) up to 6 are as follows [J: np: 0.1482, 0.0687, 0.03926, 0.0253, 0.01765; nd: 0.06358, 



0.037405, 


0.0245, 0.01724; 


nf- 


0.03524, 


0.02348, 


0.016708; ng: 


0.02237 


0.016099; 


nh: 0.015455. 


State S 






—Energy 




State 


s 




-Energy 




This work 




Literature 






This work 


Literature 



2p 

lp 
!~rp 
(yp 
3d 
id 



0.002 

0.12 

0.148 

0.01 

0.06 

0.068 

0.03 

0.039 

0.01 

0.025 

0.01 

0.017 

0.005 

0.0635 

0.005 

0.0374 



0.12300007948 
0.01747645570 
0.00009780662 
0.04561104138 
0.00447257513 
0.00030452301 
0.00503284729 
0.00010681462 
0.01040587771 
0.00012797171 
0.00467631494 
0.00027324316 
0.05056063169 
0.00005036825 
0.02626968434 
0.00000260255 



0.123000" 

0.01748858 6 ,0.01747645565 c 

0.00009780 c 

0.045611" 

0.00447257511 c ,0.004472 d 

0.005033 d 

0.010401" 

0.004651" 

0.050561" 
0.0000503682 c 
0.026270" 
0.00000260256 c 



5d 
6d 

4/ 
5/ 
6/ 

%/ 
% 
6/1 



0.005 

0.0244 

0.005 

0.0171 

0.008 

0.0352 

0.008 

0.0234 

0.008 

0.0167 

0.005 

0.0223 

0.005 

0.016 

0.005 

0.0154 



0.01505072772 
0.00004715572 
0.00899421353 
0.00006413769 
0.02330635503 
0.00002734812 
0.01216554372 
0.00004667366 
0.00624761715 
0.00000418053 
0.01503173481 
0.00004909577 
0.00896866219 
0.00005984697 
0.00895010682 
0.00003895032 



0.015051" 

0.00004715571 c 

0.008994" 



0.0000273481 c 



0.00000418053 c 
0.015032" 

0.008967" 

0.008950" 



i Ref. 



a 



b Ref. 2 



c Ref. 



■ B 



d Ref. 1 



from J-matrix calculation 26]. As seen before, in all these cases again, our energy values 
are virtually identical to this method. And whenever these are not available, the present 
work produces quite superior eigenvalues compared to the other existing values. Note that, 
for all these eigenstates, we have chosen at least one 5 value, which is very close to the S c 
value, as difficulties are encountered in these areas, with some of the methods. Near the 
critical values of 5, in general, one needs to extend r max to some larger values. For example, 
for a converged result for 4d state at 5 = 0.0374, an r max of 1100 a.u., was employed. This, 
again, is similar to a situation we came across for higher states in the previous paragraph, 
and also for some other central potentials [35 ] . 

Next in Fig. 1, variation of energy eigenvalues with respect to screening parameters are 
depicted for all the states belonging to n = 7, 8 (left) and n = 9, 10 (right) respectively, in the 
neighborhood of zero energy. Energy values increase monotonically with 5; for each n, they 
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FIG. 1: Energy eigenvalues (a.u.) of the ECSC potential for (a) n = 7, 8 and (b) n = 9, 10 levels 
respectively as a function of J in the vicinity of zero energy. 



make a distinct family and for a particular value of the quantum number n, the separation 

between states with different values of I tends to increase with an increase in S. Additionally, 

in Table III, calculated eigenvalues of all states are given at selected values of 5 (0.005 and 

0.003 for n = 8 and 10 respectively). For sake of completeness, the available S c values 

for n = 8 are mentioned in column 1 in parentheses. Only a two-parameter variational 
I I 

calculation [15j has been reported for the n = 8 state. While these are reasonable first 
estimates in absence of any other result, our GPS results are significantly better than these. 
And for n — 10, there are no results to quote for direct comparison, and it is hoped that 
these would be helpful for the purpose of future referencing. 

Now we turn to the GESC potential. Table IV reports energies for some low- and high- 
lying s states of the same for some representative b values, keeping a fixed at 1. For lower 
states, several b values are considered to understand the dependence on potential parame- 
ters. Reference results are much scarce in this case compared to the ECSC potential. The 
1/iV-expansion up to 14 terms [30] were obtained for ground and first excited state energies 
and wave functions. While these are reasonable initial estimates, clearly improved energies 
would be highly desirable. In another treatment, bound-state energies of first three s states 
have been reported within a new perturbation technique 31] . For smaller screening param- 
eters, there is, in general, a decent agreement between our result and theirs. However, the 
discrepancy starts to grow quite fast as b is increased. Finally, in Table V, eigenvalues of all 
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TABLE III: Comparison of the negative eigenvalues (a. u.) of ECSC potential for n = 8, 10 states 



at selected values of A 


. Numbers in 


the parentheses 


in column 1 denote 


5 C values 


26]. 


State 


8 


—Energy 


State 


5 


-Energy 




This work 


Literature [15] 


This work Literature 


8s(0.0100) 


0.005 


0.00314139349 


0.003134 


10s 


0.003 


0.00217587556 


8p(0.0099) 




0.00313602921 


0.003128 


lOp 




0.00217402019 


8d(0.0098) 




0.00312524768 


0.003118 


lOd 




0.00217029874 


8/(0.0096) 




0.00310894203 


0.003102 


10/ 




0.00216468974 


8g(0.0094) 




0.00308694949 


0.003080 


lO.y 




0.00215716069 


8/i(0.0092) 




0.00305904829 


0.003053 


10ft 




0.00214766774 


8i(0.0089) 




0.00302495321 


0.003020 


10/ 




0.00213615531 


8fc(0.0086) 




0.00298430926 


0.002981 


10*; 

10/ 
10m 




0.00212255542 
0.00210678692 
0.00208875461 



the I 7^ states having n < 8 are reported for the first time, for two values of b parameter. 
No results could be found for such states in the literature for comparison. This dependence 
of GESC eigenvalues on b is pictorially shown in Fig. 2 for all the states belonging to n = 2, 3 
(a), n — 4, 5 (b), n — 7, 8 (c) and n — 9, 10 (d) respectively. In all cases, energies gradually 
increase and then tend to assume a constant value. For smaller n in (a), (b), it is seen 
that, all the states belonging to a particular n form a characteristic family of the curve. 
Moreover, the states corresponding to a given n do not mix with the states with a different 
n. However, this scenario changes dramatically as we go for higher n. Thus, as we move to 
(c), appreciable complex ordering and inter-state mixing is observed for n = 6, 7 (7i mixing 
with 8s) and as we finally reach n = 9, 10 in (d), we encounter heavy mixing among the 
9i, 9k, 91 and 10s, lOp, 10<i states at around b = 0.03 — 0.06, making accurate determinate 
of these eigenvalues more and more difficult. Such complex level crossings have also been 
observed earlier for Hulthen and Yukawa potentials 35]. 



IV. CONCLUSION 

Accurate bound states of ECSC and GESC potential have been presented by means of 
a GPS method. For both these cases, zero and non-zero angular momentum states are 
calculated easily with high accuracy. The methodology is simple, efficient and, as shown, 
produces eigenvalues and eigenfunctions of comparable accuracy to those of the best available 
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TABLE IV: Comparison of the negative eigenvalues (a. 
at selected values of b. 



u.) of GESC potential for several s states 



State 


b 




-Energy 




b 




-Energy 




This work 


Literature 




This work 


Literature 


Is 


0.001 


1.99900000049 






0.005 


1.99500006211 






0.02 


1.98000390178 


1.9800039° 


,i.980oo<> 


0.05 


1.95005876574 


1.9500586" 




0.08 


1.92023217638 


1.9202305" 


,1.92023'' 


0.2 


1.80316184099 


1.8030143" 




0.1 


1.62057014563 


1.6169173" 




0.7 


1.38422179244 


1. 3477860", 1.384 6 




1 


1.19419978389 


1.0989583" 


,1.194 6 


2 


0.82070036307 






3 


0.66846103237 






5 


0.56680152293 






10 


0.51787565892 






20 


0.50467744871 




2s 


0.001 


0.49900000697 






0.005 


0.49500085803 






0.02 


0.48005182598 


0.4800516" 


,o.480oo<> 


0.1 


0.40487183925 


0.4043555",0.4048 6 




0.3 


0.27382625160 


0.2431595" 


,0.274 6 


0.5 


0.21294420503 






1 


0.17216986942 






5 


0.13273550025 






10 


0.12719012732 






15 


0.12601285391 




3s 


0.001 


0.22122225644 






0.02 


0.20245702303 


0.2024526" 




0.06 


0.16695326252 


0.1662097" 




0.2 


0.09999056593 






0.5 


0.07573455915 






1 


0.06809011302 




4s 


0.001 


0.12400010645 






0.05 


0.08216008138 






0.1 


0.05956143168 






0.5 


0.03910015933 




5s 


0.001 


0.07900025672 






0.2 


0.02693147558 




6s 


0.001 


0.05455608155 






0.1 


0.01998364734 




9s 


0.001 


0.02369391160 






0.1 


0.00776025524 





^Ref. 



a 



b Rcf. 3i 



methods found in the literature. All the 55 states lying with n < 10 for the former, and 
36 states with n < 8 for the latter, are calculated up to eleven significant figures covering 
wide ranges of interaction. For the former, our results are superior to all the existing 
methods except that of the J-matrix formalism, while for the latter potential, our results 
surpass the accuracy of all existing methods. A detailed analysis of the variation of energies 
with respect to potential parameters show quite different trends for these two potentials. 
For higher n, complex level crossing and inter-state mixing has been observed for the GESC 
potential. Special attention was paid for the high-lying states and regions of strong screening 
parameters. Many states are presented here for the first time. This offers a simple reliable 
method for the accurate calculation of these and other potentials in quantum mechanics. 
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TABLE V: Energies (a. u.) of the GESC potential for several / 7^ states at selected values of b. 



State 


6 


—Energy 


b 


—Energy 


b 


-Energy 


2p 


0.001 


0.49900000498 


0.2 


0.32215735500 


1 


0.13471500886 


3p 


0.001 


0.22122225198 


0.05 


0.17482725023 


0.5 


0.06796045435 


3d 




0.22122224308 




0.17414406972 




0.05712005483 


■\p 


0.001 


0.12400009861 


0.05 


0.08176332315 


0.3 


0.03976138004 


Ad 




0.12400008291 




0.08094164047 




0.03494876135 


4/ 




0.12400005934 




0.07963728134 




0.03166956721 


5p 


0.001 


0.07900024460 


0.05 


0.04265374603 


0.2 


0.02588418263 


5<i 




0.07900022035 




0.04187154574 




0.02383314014 


5/ 




0.07900018393 




0.04063035100 




0.02127952212 


% 




0.07900013528 




0.03883963267 




0.02015250137 


ftp 


0.001 


0.05455606434 


0.05 


0.02479616384 


0.2 


0.01712467704 


6d 




0.05455602989 




0.02418563145 




0.01601163848 


6/ 




0.05455597816 




0.02322582411 




0.01465426827 


6c/ 




0.05455590906 




0.02185943895 




0.01401444073 


6/1 




0.05455582250 




0.01999675997 




0.01389838294 


Ip 


0.001 


0.03981726562 


0.02 


0.02456106851 


0.1 


0.01354669295 


Id 




0.03981721949 




0.02443236181 




0.01300727690 


7/ 




0.03981715020 




0.02423479689 




0.01221827987 


7g 




0.03981705766 




0.02396283009 




0.01127101023 


Th 




0.03981694172 




0.02360879174 




0.01052800497 


1% 




0.03981680223 




0.02316253952 




0.01025534957 


Hp 


0.001 


0.03025158916 


0.01 


0.02226047376 


0.1 


0.00996391701 


8d 




0.03025153001 




0.02223011222 




0.00961764441 


8/ 




0.03025144118 




0.02218403767 




0.00911555429 


89 




0.03025132253 




0.02212160821 




0.00852094708 


8/1 




0.03025117390 




0.02204196197 




0.00805230527 


Si 




0.03025099505 




0.02194401054 




0.00786434715 


8fc 




0.03025078574 




0.02182642900 




0.00781852403 
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FIG. 2: Energy eigenvalues (a.u.) of the GESC potential for (a) n = 2, 3 (b) n = 4, 5 (c) n = 7, 8 
(d) n = 9, 10 levels, respectively, as a function of b in the vicinity of zero energy. 
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